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o 

. We study the thermodynamics of the linear sigma model with constituent quarks beyond the mean- 

fSJ ' field approximation. By integrating out the quark degrees of freedom we derive an effective action 

for the meson fields which is then linearized around the ground state including field fiuctuations. We 
propose a new method for performing exact averaging of complicated functions over the meson field 
fiuctuations. Both thermal and zero-point fiuctuations are considered. The chiral condensate and the 
effective meson masses are determined self-consistently in a rigorous thermodynamic framework. At 
zero chemical potential the model predicts a chiral crossover transition which separates two distinct 
regimes: heavy quarks and light pions at low temperatures, but light quarks and heavy mesons at 
high temperatures. The crossover becomes a first order phase transition if the vacuum pion mass is 
reduced from its physical value. 
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o 

CN ■ I. INTRODUCTION 

o ■ 

(j~>) , It is commonly accepted that Quantum Chromo-Dynamics (QCD) is the true theory of strong interactions. There- 
fore, in principle, it should describe all phases of strongly interacting matter at all densities and temperatures. In 
practice QCD can be exactly solved only in some limiting cases: firstly, at very high densities and temperatures 
when the property of asymptotic freedom allows a perturbative expansion to be used, and secondly, at zero density 
O and high enough temperatures when it can be discretized on an Euclidean lattice. Although some improvements at 
^ finite chemical potential have been achieved recently 0, the direct application of QCD at high baryon density, or 
^ for real-time processes, is still quite problematic at present. This makes it necessary to build effective models, which 
^ , respect only general properties of QCD, such as chiral symmetry, and operate with effective degrees of freedom, such 

' mesons and constituent quarks. 
rS [ The linear sigma model (LctM) |3_is one of the most popular models of this kind which has been studied already 
■ for several decades (see e.g. Refs. 0, Q). It incorporates correctly the phenomenology of low-energy strong inter- 
actions, including chiral symmetry. At low density and temperature the matter is assumed to be in a phase where 
chiral symmetry is spontaneously broken. A phase transition that restores chiral symmetry was predicted at high 
temperatures 5] or baryon densities According to Ref. QCD with two massless flavors belongs to the same 
universality class as the 0(4) LaM and therefore the phase transition is of the same order. In the case of nonzero 
quark masses chiral symmetry is explicitly broken and universality arguments cannot be applied. Then the character 
of the chiral transition depends on the detailed pattern of the symmetry breaking Q . In recent years there have been 
many attempts to use different versions of the LcM to model QCD phase transitions at finite temperature 0, H, j 
baryon density 0, ^| , and for non-equilibrium conditions 0, 0, IT^ Hsl IT6l | . 

Despite many studies the status of the LctM remains somewhat controversial. In most applications of this model 
only mesonic degrees of freedom are included explicitly 0, 0, ^3 ■ Then the model shows a chiral phase transition, but 
the high-temperature phase is very different from what is expected for QCD, since there are no free quarks but only 
heavy mesons. On the other hand, if quark degrees of freedom are included from the beginning while fluctuations in 
the mesonic fields are ignored [TH . the model also predicts a chiral transition at about the same temperature. Now, 
however, the low-temperature phase has the wrong structure since it is dominated by constituent quarks and, in the 
mean field approximation, pions play no role. Physically, of course, pions are expected to be the most relevant degrees 
of freedom at low temperatures. There have been several attempts to improve matters by including both constituent 
quarks and pions (see e.g. Ref. ^3)' but the calculations have been limited to the lowest order loop expansion. 
A more satisfactory approach has been pursued recently, namely the (approximate) solution of the renormalization 
group flow equations 17, 18] which include effects due to field fluctuations. This appears to be particularly valuable 
in the neighborhood of critical points. 
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In the present paper we study a different, and possibly more transparent, approach to including field fiuctuations. 
We deal with the full LcrM including both constituent quarks and mesonic excitations and our goal is to proceed as far 
as possible without invoking any kind of mean-field approximation or perturbative expansion. We shall demonstrate 
below that one can indeed develop a practical computational scheme where the field fluctuations are incorporated in 
the thermodynamic potential to all orders in a self-consistent way. This formalism corresponds to summing up the 
infinite series of daisy and superdaisy diagrams. Of course, this is possible only within the Hartree approximation 
where the exchange diagrams are disregarded. The resulting physical picture appears to be close to QCD-based 
expectations, e.g. Ref. [13 ■ Namely, the model exhibits a smooth chiral crossover transition at temperatures of 150- 
200 MeV from a low-temperature phase made of heavy constituent quarks and light pions to a high-temperature phase 
composed of light quarks and heavy mesonic excitations. We regard these results as quite satisfactory for modeling 
QCD, even though the model has no confinement. For effective theories that incorporate both chiral symmetry and 
confinement see Ref. pol |. 

The paper is organized as follows. In Section II a general LcrM is formulated and its partition function is represented 
within the path integral formalism. We integrate out the quark degrees of freedom and approximately reduce the 
problem to a purely mesonic theory with a very non-linear effective potential. This potential is then linearized around 
the correct ground state including equilibrium field fluctuations to all orders. Finally, the mesonic contribution to 
the thermodynamic potential is calculated in closed form. In Section III we describe a general method for evaluating 
averages of complicated functions over the fluctuations of the fields. This permits easy consistency checks and allows 
the equations to be put into a simple form. Results obtained by including or excluding zero-point contributions, in 
addition to the thermal fiuctuations, are presented in Section IV. Here we also assess simplified approximations to 
the complete results. Our conclusions are presented in Section V. The relationship of the present approach to earlier 
work |2^|23| is discussed in Appendix A, and the proof of a useful identity is given in Appendix B. 



II. THEORY 



A. LcrM Lagrangian and Partition Function 



In this paper we employ a standard version of the LcrM model with SU(2)l x SU{2)ji symmetry. The corresponding 
Lagrangian for quarks interacting with the a and tt meson fields is written as 



where the quark Lagrangian is 



the meson kinetic energy is 



and the meson potential is given by 



Cqq + CKm-U{a,-K) , (1) 

^qq = q [ilfj^d'^ ~ 9(0- + ilbT ■ tt)] q , (2) 

CKm = \ {d^Td'^a + d^TV ■ a^TT) , (3) 

Uia,7r) = \X{a' + n'-(:)' -Ha. (4) 



Here we have included an explicit chiral symmetry breaking term of the conventional type. Ha. As is well known, the 
choice of the symmetry breaking term is not unique. The parameters of the model, g, A, C and H, will be specified 
later. The partition function of the grand canonical ensemble can be written as a functional integral over the quark 
and meson fields [23 



Z = Trexp |-/3(i7 - ^iV)| = J VqVqVaVn cxp I J dr J (fx{C 

1 V 



mi\) ) , (5) 



where TV is the quark number operator, (3 = T ^ is the inverse temperature, /i is the quark chemical potential, V is 
the volume of the system and t = it denotes the imaginary time. 
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B. Integrating Out the Quarks 

First we integrate out the quark degrees of freedom. The relevant part of the partition function is 

/3 



Zqq 



JvqVqexp I J dr J cPxqDq \ , (6) 



. V 



where 



d 

D = + 17 • V - .g(cr + i75T • tt) + /i7" . (7) 

Formally one can integrate out the quark fields with the result 

InZ^, = j dr j d^xlndeti) . (8) 

V 

The analysis of such an expression has been discussed in detail by Fraser p3 |. Accordingly, one proceeds by moving 
the operators to the left and the fields which depend on r and x to the right. This generates a series of commutators 
which involve the derivatives of the meson fields. Since our interest is in low energy properties we will discard 
these commutator terms, tacitly assuming that the meson mode amplitudes vary slowly in position and time. Then, 
evaluating the operator 13 in a frequency-three-momentum representation we get 

Indet ^ = ^ 5] {In [/3^(^^ + {E - ^,f] + In [(3\u^l + {E + ^^f] } , (9) 

pn 

where the Matsubara frequency LOn — (2n -I- l)7rT, the quark degeneracy Vq is 12 for the two flavors employed here, 
and E^ — -\- m? . The quark effective mass is given by 

^g\a^ . (10) 

After performing the summation over n we get In Zqq — ~ J^^ dr Jy (Pxnqq{m), where the quark-antiquark thermo- 
dynamic potential density is expressed as 



nqq{m) - - ^ y dpp^ {PE + In [l + e-/3(S-p) 



In 



l + e" 



(11) 



Note that this differs from the standard result in that the mass depends on the meson fields according to Eq. I|l()|l . 
Thus the partition function lO can be written as 

Z = J VaVnexp I J dr J cPxCm ) , (12) 

I V 

where 

Cm^ CK,n~U{cr,T^) (13) 

is the effective meson Lagrangian. Here the effective potential, including the contribution of quarks and antiquarks 
with effective mass to, is 

U{a, tt) = U{a, n) + nqq{m) . (14) 

The Lagrangian H13|l now contains only meson fields and it constitutes an effective meson theory with the very 
nonlinear interaction potential H14|) . 
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C. Linearization of the Mesonic Action 



Let us consider the equations of motion for meson fields which follow from the effective meson Lagrangian H13() : 

dU 

d^d^a + ^ - , (15) 



d^d^TT + — = . (16) 



dU_ 

(57? 

We average these equations over the meson field fluctuations. Since the a field is expected to develop a nonvanishing 
expectation value u, we decompose it as cr = + A where A is the fluctuating part. By definition the average of A is 
zero and this is also true for any odd power of A, i.e. (A") = for odd n. A similar remark applies to the pion field. 
Here and below angle brackets indicate averaging over the field fluctuations. A practical scheme for evaluating such 
averages is discussed in the following section. 

As the thermal average of an odd number of fluctuating fields is zero, only the terms in U with odd powers of A 
will contribute to Eq. (|15|l . yielding 

Since the pion field always occurs as tt^, a single derivative with respect to a component tt^ will always yield an odd 
number of fluctuating fields and thus the thermal average automatically vanishes. Notice that H17|) includes all powers 
of the fluctuating fields, so the equation of motion contains fluctuations of the flelds to all orders. This makes our 
approach different from numerous previous attempts to include fluctuations where an expansion is made around the 
mean field ground state. The second derivatives of the effective potential pick out the even powers of the fluctuating 
fields in U and these we identify with the meson masses: 

We now linearize the complicated effective mesonic potential by setting 

U{v + A,7r) ^ (U{v + A, tt)) + im^ (A^ - (A^)) + im^ (tt^ - (tt^)) . (19) 

Obviously, this becomes an identity if we average over the field fluctuations on both sides. The terms on the right 
containing the average quantities should be included directly in the total thermodynamic potential density. The 
remaining terms containing A^ and tt^ are combined with the kinetic energy to give the mesonic partition function: 

Z,„ = / PApTrexp <( / dr / d^x [Ck^ - ^m^A^ - im^vr^] ) . (20) 

Following standard steps one arrives at the thermodynamic potential density associated with the meson field 
fiuctuations: 

pV 




T 

2^ 



J dpp" {\(3E, + In (1 - e-^^')] , (22) 
" S / '^^^^ ^^'^^^ + " ^ ' ^^^^ 



where E„ ~ \/p'^ + "m-a — + mj. Two consistency relations for the meson masses follow directly from 

Eq. (1201): 

(A")^2|| , (,^).2|2|. (24, 

Finally, we can write the total thermodynamic potential density as 

n = {U{v + A,7r)) + (ri,-,(m)) - \ml{^^) - h^mHn^) + a„(m,,m,) . (25) 
The subtraction of the third and fourth terms on the right is necessary to avoid double counting psf. 
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III. EVALUATION OF FIELD FLUCTUATIONS 



A. Averaging Procedure 

Consider now the average of a complicated function 0{v + A, tt^) over the fluctuating fields A and tt. The averaging 
can be carried out by generalizing the technique introduced in 26]. First, we expand the function in a Taylor series 
about the point {v,0), then take the average term by term 

/ A fe 2n \ 

(O(. + A,.r^))=;^O('=-)(z;,0)(^^) , (26) 

k,n \ ■ ■ / 

where 

nfe an 

(a, b) = — — 0(a, b) . (27) 

The next step is to reduce a general vertex (A'^tt^") to powers of (A^) and (tt^). The necessary counting factors for 
joining the meson fields at this vertex in all possible ways are: (A*"') = (k — 1)!!(A^)'^/^ for k even, and zero for 
k odd. For pions, all species are equivalent so (nf) ~ (7r|) — (7r|) = |(7r^) and therefore (tt^") = (2n + l)!!(i7r^)". 
After substituting these factors in the series H26|l we notice that the resulting averaging is equivalent to performing 
integrations with the following Gaussian weighting functions: 

P4z) = (27r(A2))-V2exp(^_-^^ , (28) 
P^y) - \fl ( A)^ . (29) 



Then, resumming the Taylor series, we obtain 

{0{v + A,7v^))= I dzP^z) I dyy'PAy)0{v + z,y') . (30) 



This is a general result for any analytic function O. Note that {0{v + A, tt^)) 0{v, 0) in the limit when (A^) and 
(tt^) 0. The correspondence between Eq. l|5n|l and the approximate expressions used in Refs. [^H^l is discussed 
in the Appendix. We also need the derivative of Eq. (|Sn|) with respect to some variable a. After two integrations by 
parts one obtains 

-(0(. + A,7r )) = ^^ g-^ /+2^\ 5A^ /+2^\ d^f / " ^^'^ 

Using this equation one can easily check that the derivative of the total thermodynamic potential density with 
respect to v is 




= . (32) 

This vanishes due to the equation of motion H17|l and the mass definitions (|18|l . Thus we have the non-trivial result 
that f2 is a minimum with respect to variations in the scalar condensate v, as is required for sensible thermodynamics. 
This analysis amounts to the formal justification of the approach of Refs. pH l2a | where the thermal averages were 
evaluated approximately in series form, and therefore thermodynamic consistency was obtained only approximately. 



B. Zero-point Fluctuations 



The general formulae 122|l and H23|l include both zero-point and thermal fiuctuations. The former are divergent 
and thus require a proper renormalization procedure. The non-trivial issue of renormalization in self-consistent 
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approximation schemes has been discussed for the LctM model by several authors |2 
for the regularization in four dimensions, namely 



^|29|. We adopt their result 



di7vi 



1 

3 
4^ 
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ml 




^ 167r2 


A2 
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ml In 


ml 




167r2 


A2 

m 



A 



(33) 



Note that the +A^ term may be removed by redefinition of the constant ( in the sigma model potential The 
zero-point contributions to then follow by integration (see Eqs. (|35|l and H36|l below). 
For the quark case it is natural to adopt a similar form, setting 



c?m2 



27r2 



dp 



87r2 



'In 



A2 

9 



A2 



(34) 



As is reasonable in a phenomenological model we have allowed the renormalization scale Aq in the fermion sector to 
differ from the scale A^ in the meson sector to allow for the different physics in the two sectors. Again, integration 
gives the result for riqq (see Eq. below). 



C. Practical Evaluation 



The thermodynamic potential densities including both zero-point and thermal fluctuations of the sigma and pion 
fields are, respectively: 



1 


In 


ml 


647r2 


A2 


3 


In 


ml 


647r2 


A2 



4-5/dp/ln(f-e-^- 



The squared fluctuations of the meson fields are found by using Eq. H24|l : 



1 


ml In 


ml 


167r2 


A2 

m 


3 


ml In 


ml 


167r2 


A2 



1 f p2 

/ dp—nB{E„) , 
J E^ 



^l^-ml 



(35) 
(36) 

(37) 
(38) 



where the Bose occupation number is nsix) — (exp f3x — 1) ^. The thermodynamic potential density associated with 
quarks is: 



327r2 



1 + e-'^^^-^) 



In 



(39) 



where E = \/p'^ + m^ . 

Since the total thermodynamic potential (j25|l is a minimum with respect to variations in the scalar condensate v, 
the thermodynamic quantities can be found in standard fashion. The pressure P = — O and the energy density is 



£ ^ ({/(7; + A,7r))-X(^')-X(^' 



1 



3 


m^ In — 

m 


647r2 




r Af 


327r2 



647r2 

2 A 2 ^,- A 2 ™2 ^ 



m^ln^ + i(3m2-A^)(A^^-m^) 



^ j dpp^ [E^nsiE,) + SE^nsiE^)] 



i_ In — + i(3m2 - A2)(A2 - m^)^ + ^ / '^PP' (i?["F(£;, ^i) + ^^(ii;, -^)]) , (40) 
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where the Fermi-Dirac occupation number is np{x,y) = [exp/3(a; — y) + 1]^^. Here we do not expUcitly show a 
constant which must be subtracted from the pressure and added to the energy density in order to render P and £ 
zero in the vacuum. The quark number density is 



dvp" {nF{E,^i)-nF{E,-^i)) . (41) 



The entropy density can then be obtained from the standard thermodynamic relation S — f3{£ + P — /in). 
For the equation of motion (|17|l we need the first derivative of flqq : 



99 



dA 



.9' (aA(m)) , (42) 



where the function A(m) is defined by 
Mm) = 2ff = 



2 

2 1 A 2 2 

m In —rr + A„ — m 

g 



/2 
dp^[nF{E,^i)+nF{E,-^i)]. (43) 



Within our approach the quark condensate can be easily expressed as 



g \ da 

For the meson masses in Eq. H18(l we need second derivatives 



{m)--(^)=9{<jAim)) . (44) 



d'i^,A ^^2/^(^^)^25V^^\ = ^(A(« + A)^(m)) , (45) 



where for the latter equalities we have used Eq. H30() . Explicitly 

^ — ^ In - - ^ y d,-[nFiE^ ,) + UFiE^ -,)] • (47) 

When deriving Eqs. (|43|l . I|45|l and H46|l we have used the fact that Vlqq{m) and Aim) are even functions of m. For 
simplicity we have suppressed the dependence of these functions upon T and /i. 

Finally, we need to average the bare potential [/(v + A, tt) and its derivatives. Using the expression H3(J|I . or more 
elementary means, one obtains 

(t/(z; + A,7r)) - (iA(a2 + 7r2-C)'-i/a)-iA{3(«2 + (A2»2 + („2 + (^2))2 

+ |(7r2)2 + 2(A2)(7r2) - 2C(v2 + (A^) + (tt^) - !() - Sw^} - iJv . (48) 

The first derivative of U needed in the equation of motion (|17() is 

C^^V/'"^ ) = Mv^ + 3(A^) + (.^) -C)-H. (49) 

The second derivatives of U needed for the meson masses in Eq. (|18ll are 

'a2i7(w + A,7r 



5A2 

92[/(t; + A,7r) 



= \{?,v' + 3(A2) + (tt^) - C) , 

= X{v^ + {A^) + U^^)-C). (50) 



The equation of motion p7|) involves the sum of Eqs. (|42|) and (|49|) . It is interesting to note that the quark 
contribution of Eq. H42|) vanishes if w = 0. This is because {aA{m)) becomes {AA{g^/S^^\^^z^)) which is an odd 
function of A. In the absence of explicit symmetry breaking, H — 0, Eq. I|49|l permits a solution v — which is the 
true solution above the critical point. 
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TABLE I: Model Parameters 



parameter 


ZPT 


NOZPT 


A 


7.114 


27.58 


C (MeV^) 


-8.733 X 10'' 


7.847x10^ 


H (MeV^) 


1.760x10*^ 


1.760x10^^ 


g 


2.844 


3.387 


Am (MeV) 


453.4 




Ag (MeV) 


951.8 





The equation of motion and the equations for the meson masses have to be solved self-consistently. The integrals in 
Eq. (|30|l were evaluated numerically using 32-point Gaussian integration |30l |. The necessary thermodynamic integrals 
were obtained by using the numerical approximation scheme of Ref. |3]| . 

Assembling the contributions to the meson masses in Eq. H18|) and using the equation of motion H17|) for the case 
?; ^ 0, we find 

ml ^ 2Xv' + ^ ' l) + + f ' (51) 

ml = 2A {{nl) - (A^)) + 5^ ( " ^) Mm)) + f • (52) 

It is interesting to consider the high temperature limit when (A^) and (Trf) become equal. In this limit the first term 
on the right in eq. H52|l vanishes, but so does the second quark term too (see Appendix B). Therefore, we regain the 
usual high temperature result, ~ H/v. Since in this region v is very small, the quark contribution in (|51|l also 
vanishes approximately and the first term in the equation can be neglected, returning ~ TOtt as expected. 

IV. RESULTS 
A. Choice of Parameters 

We need to choose the parameters A, Cj g, Am and Ag in the case where zero-point fluctuations are included 
(labeled ZPT). We also present calculations where zero-point fluctuations are excluded (labeled NOZPT) for two 
reasons. Firstly they have often been excluded in the literature and secondly it is important to assess their influence 
in the ZPT case. Note that the parameters A^ and A^ are irrelevant to the NOZPT analysis. 

Since chiral symmetry is spontaneously broken in the vacuum, the axial current requires that Uvac = /ttj where the 
pion decay constant is f-^ = 92.4 MeV. Goldstone's theorem states that in the absence of explicit symmetry breaking, 
H = 0, the pion, as a Goldstone boson, should have zero mass in the phase with spontaneously broken symmetry. 
Then Eq. I|52|l gives in the vacuum 

2A {{nf) (A^)) + g' ^ ^) A{m)) = . (53) 

It is natural to assume that this equation is also valid in the case ^ 0, or, stated differently, to take H = tti^/tt. 
The solution of Eq. (|53|l at T = is (A^) — (tt^) since, as shown in Appendix B, the quark contribution vanishes in 
this limit. This same condition was used by Lenaghan and Rischke [231 in the pure meson sector. For the NOZPT 
case (A^) — (tt?) = in the vacuum so the equation is automatically satisfied. In the ZPT case Eq. H53|l requires Am 
to be intermediate between the pion and sigma masses. Choosing vacuum masses = 700 MeV and — 138 MeV 
gives Am = 453 MeV, as listed in Tabled values of similar order have been employed in Refs. [28lf29ll3^ . This should 
not be directly viewed as a cut-off parameter. If the m\ term in Eq. 1)35(1 or the term in Eq. 1(37(1 is compared 
with the result of using a four-dimensional cut-off to evaluate the integral, the appropriate cut-off parameter is more 
reasonable, namely A^e^^^ — 960 MeV or AmC^^^ — 748 MeV respectively. The value of the coupling constant g is 
fixed by the requirement that the vacuum quark mass be approximately a third of the nucleon mass. We can define 
the mass by M„ = {(to")}1/". It is not a 'priori obvious which power of the mass should be averaged. We examine 
the two lowest moments with n = 1 and n = 2 which are 

Ml = (m) =5(v/^^T^) , (54) 

Af2 = (m2)5 ^ g(y^ + (A2) + (Tr^))^ . (55) 
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FIG. 1: The a mean field, i;, and the root mean square fiuctuations as a function of temperature (a) without and (b) with 
zero-point fluctuations. 

As we shall see below, the two definitions give masses which are closely similar, so it is of little consequence which 
choice is made. We actually used the second definition, setting M2 ~ 939/3 MeV. To determine A and Kq in the ZPT 
case we require the vacuum mo- to be 700 MeV in Eq. I|51ll and constrain the vacuum quark condensate in Eq. H44|l . 
For the latter we choose the value (917) vac = {uu + dd) = -2 x (225 MeV)^ HIQ. In the NOZPT case the vacuum 
quark condensate is zero and only the a mass equation is required. The value of the scale Ag in Table^is about twice 
that in the meson sector, as seems intuitively reasonable. We remark that if Ag is chosen to be equal to A^, and the 
quark condensate is predicted, rather than fitted, it turns out to be more than an order of magnitude smaller than 
the physical value used above. Furthermore, even when these renormalization scales differ, it is not possible to carry 
out the fit outlined above with rriu = 600 MeV. This is why we have chosen a somewhat larger value of 700 MeV for 
this poorly known mass. Thus, the vacuum value of the quark condensate is a strong constraint. It is also responsible 
for the quite different values of ( listed in columns two and three of Table These were obtained from the equation 
of motion (|17|l . 

Finally, we mention that we have examined results with a number of different parameter sets and qualitatively they 
are all very similar. Therefore the parameters of Table |l] will provide a representative set of results for this model. 

We should further point out that while we require Goldstone's theorem to be satisfied at T = 0, this docs not 
guarantee that it will be fulfilled in general. This is because for T, /i 7^ the r.h.s. of Eq. (|52|l does not vanish even 
when the explicit symmetry-breaking term H is set to zero. The nonvanishing terms can be cancelled by including 
an additional (exchange) diagram [Hill which is missing in our treatment. However, we focus on the more realistic 
case where the pion has its physical vacuum mass, rriTr = 138 MeV. Since we find that the contribution of these 
nonvanishing terms is relatively small, we shall keep the form 1)52(1 for the pion mass so as to preserve the structure 
of the theory. 



In this paper we consider the case where the quark chemical potential, /i, is zero, corresponding to a net quark 
density, n, of zero. The behavior of the sigma mean field v is shown for the NOZPT case in Fig. ^a) and for the ZPT 



B. Full Model 
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case in Fig. ^b). The value of v decreases smoothly to zero as the temperature increases so that chiral symmetry is 
approximately restored at high temperatures. Thus, there is no sharp phase transition, but rather a crossover. The 
crossover temperature, determined from the maximum value of \dv/dT\, is 150 MeV and 198 MeV in the NOZPT 
and ZPT cases, respectively. We remark that if quark contributions are omitted this temperature is about 70 MeV 
higher in accordance with Ref. |23|, and in that case v falls off less rapidly with temperature. Lattice calculations 
suggest a smooth crossover, as we have found, and extrapolation to the chiral limit yields a phase transiiton with a 
critical temperature of 173 ±8 MeV (19|]. The phase transition is thought to be of second order, although first order is 
a possibility Our model is unlikely to be accurate enough to distinguish the order, but it is of interest to examine 
the qualitative predictions as we approach the chiral limit. We do this in the simplest fashion by multiplying the 
physical vacuum pion mass by a factor < 1 (the ZPT parameters are refitted with the remaining quantities unchanged, 
as specified in Subsec. IVA). We find that when the pion mass reduction factor is | the smooth crossover becomes 
a first order phase transition with a critical temperature of 197 MeV. This temperature is little changed by further 
reduction in the pion mass. For example with a reduction factor of j the temperature is 194 MeV. This appears to 
be an improvement on the prediction of the LtrM without quarks [3^ since it agrees with the lattice finding of a weak 
dependence of the critical temperature on the quark mass UJ. In comparison the renormalization group calculations 
\IT\ show a somewhat stronger dependence on the pion mass and further a phase transition (second order) is obtained 
only in the chiral limit. The transition temperature of 100-130 McV found in Ref. ^3 is substantially below the 
values given above, no doubt due to the rather small value of 430 MeV obtained for the a mass. 

Figuren]also shows the mean square fluctuations, (A^)^ and (7rf)2. We required them to be equal at T = (see 
previous subsection), and they are equal again at high temperatures when v becomes small. As can be seen from Fig. 
^a) the thermal contribution in the intermediate region is larger for the pion since it has the smaller mass. We do 
not see the significant enhancement of the field fluctuations in the crossover region that would be expected were there 
to be a second order phase transition. The reason is that large thermal fluctuations are developed only at <SC T, 
which never holds in our calculations: in fact Fig. |21 shows that TOct docs not drop below 300 MeV. The shallow dip 
in (A^)2 in the 100-200 MeV temperature range which is seen in Fig. Q^b) arises from the decrease in mo-. This 
causes the zero-point contribution to the fluctuations in Eq. H37|l to decrease, becoming zero at mo- = A^. At higher 
temperatures Too- rises with increasing T, causing the zero-point contribution to increase strongly and dominate over 
the thermal contribution for T > 300 MeV. This is clearly seen by comparing Figs. H^a) and (b). 
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FIG. 3: The thermodynamic quantities SP/T"*, E jT^ and as a function of temperature (a) without and (b) with zero-point 

fluctuations. 



The behavior of the meson and quark masses is illustrated in Fig. |21 For the meson masses the NOZPT and 
ZPT calculations give qualitatively similar results. As expected, the pion mass is an increasing function of T while 
the sigma mass develops a minimum in the transition region. Here the pion and sigma masses become nearly equal 
signaling the approximate restoration of chiral symmetry. At high temperatures, T > 250 MeV, the meson masses 
become large and grow according to H/v. For the quark mass it is seen from Fig. [51 that the definitions H54|l and 
()55|l for Ml (dotted curve) and M2 (short-dashed curve) give closely similar results. At high temperatures the mass is 
small in comparison to the temperature in the NOZPT case (Fig. |2ta)), as it is in the renormalization group approach 
[iTj . In contrast in the ZPT case the increase in the fluctuations seen in Fig. ^ causes an increase in the quark mass, 
such that it is of the same order as the temperature. Caldas et al. js^ report a similar result in their approach to 
including fluctuations. In this high-T regime the meson masses in Fig. EJb) are approximately twice the quark mass, 
as would be expected for loosely bound states of quarks and antiquarks. It should be emphasized that this behavior 
is qualitatively different from the mean- field approximation where the constituent quark mass is defined to be gv^ so 
that it vanishes at high temperatures. 

The thermodynamic quantities, energy density, pressure and entropy density, as functions of temperature are 
presented in Fig. |31 At high temperatures these are dominated by quarks. For a massless quark gas of degeneracy 
12, 2,P/T^ = f/T" = 6.91 and S/T^ = 9.21, and in the NOZPT case in Fig. O^a) these values are achieved at 
high temperatures. There is only a small contribution from mesons, since they are heavy, indicating that the system 
has effectively become a massless quark gas. The same cannot be said for the ZPT case in Fig. Efb), where the 
asymptotic values are smaller than the massless limit. This is due the quark mass remaining comparable to the 
temperature. The quark condensate is displayed in Fig. 0]for the ZPT case since the NOZPT case has only a small 
thermal contribution. The condensate is dominated by the zero-point contribution. Even though one might expect 
the thermal part to increase with T, the condensate (|44|l is approximately gv{A(m))^ so the decrease in v is sufficient 
to bring the thermal contribution to zero at high temperature. Thus {qq) starts at T = with the chosen empirical 
value and rapidly becomes negligible when v becomes small. This is the physically expected behavior. 
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FIG. 4: Thermal, zero-point and total values of the quark condensate as a function of temperature in the ZPT case. 

C. Approximations 

The expressions above which involve Qqq all require the calculation of two integrals over the Gaussian weighting 
functions as well as a momentum integration. Thus it is worthwhile to examine simpler procedures. The most natural 
simplification is to thermally average the quark mass and then insert it in Eqs. H42() . J45p and (|46|l . We first consider 
the Ml mass definition in Eq. H54|l which was the approximation adopted in Ref. ■ Replacing m with Mi, the 
thermodynamic potential density is approximated {^lqq{m)) flqq{Mi). In the equation of motion 

and for the meson masses 

The M2 case defined in Eq. H55|l is even simpler. Here the thermodynamic potential density is approximated 
(rigq(m)) — > nqq{M2), and in the equation of motion 

(^) ^ s'{'r)AiM,) = g\A{M,) . (58) 

The second derivatives are g'^A{M2) for both sigma and pion. In this case the quark contributions to the meson 
masses in (I51|) and H52|) vanish. 

Note that ambiguity in the evaluation of Eqs. (|42|l . (|45|l and l|46|l when using a thermally averaged quark mass is 
resolved by using Eq. H32|l to require that dfl/dv be identically zero. Then the thermodynamic potential density is a 
minimum with respect to variations in the scalar condensate v, as it should be. 

In Figs. OandElwe compare the full ZPT results with those obtained by using the Mi and M2 approximations for 
various quantities of interest. The parameters in the left column of Tableware used throughout. At temperatures 
above 220 MeV both approximations are able to reproduce the exact results reasonably well with, perhaps, a slight 
preference for the Afi case. Below that temperature there are larger deviations between the exact and the approximate 
results. The crossover occurs at a 15 MeV higher temperature in the Mi approximation. At lower temperatures 
the Ml approximation gives a reasonable, although not extremely accurate account of the exact results. The M2 
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approximation, on the other hand, shows marked deviations for the sigma mean field and mass, v and m^r, and the 
quark condensate. Of course, this can be much improved by refitting the parameters at T = 0. However, then one finds 
that there are no physical solutions at temperatures just above 300 MeV. We conclude that the Mi approximation is to 
be preferred and that it describes the trends of the exact results quite well, being most accurate at high temperature. 



V. CONCLUSIONS 



Our main goal in this paper was to understand the role played by field fiuctuations in effective chiral models. We 
have shown how to average field functions of arbitrary complexity over the field fiuctuations. This allowed a set of 
self-consistent equations for the average field and the masses to be formulated which led to consistent thermodynamics. 
We have applied this approach to the linear sigma model, including both mesonic and quark degrees of freedom. The 
quark degrees of freedom were integrated out and the effective action was linearized according to Eq. H19|l . Thus 
we described the a and tt mesons as quasiparticles and their properties were properly taken into account in the 
thermodynamic potential. 

We have considered two versions of the model: ZPT where both zero-point and thermal fluctuations are included, 
and NOZPT where only thermal fluctuations are present. The former was able to describe the quark condensate and 
the vacuum value was chosen according to QCD-based estimates. To accommodate this value a sigma mass of 700 
MeV was employed which, while within the broad range set by the Particle Data Group , is 100 MeV larger than 
the traditional value. We also needed to employ separate renormalization scales for the mesonic and quark zero-point 
contributions of 450 MeV and 950 MeV, respectively. For both versions of the model we required the constituent 
quark mass in vacuum to be 1/3 of the free nucleon mass. 

Numerical results were presented only for the case of zero net quark density, i.e. /i = 0. The calculations revealed 
an interesting and consistent picture. With increasing temperature we saw a gradual decrease of the condensate and 
an increase of the sigma and pion fluctuations. The model gave a crossover type of chiral transition at a temperature 
of about 198 MeV (150 MeV) in the ZPT (NOZPT) case. In the transition region the rms field fiuctuations became 
comparable in magnitude to the condensate. On the other hand we did not observe the strong increase in the sigma 
field fluctuations which would appear in the vicinity of a true critical point. The restoration of chiral symmetry was 
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FIG. 6: Comparison of the full results with the Mi and M2 approximations: (a) thermodynamic quantities "iP jT^ and E jT^ 
and (b) quark condensate. 

seen in the behavior of the sigma and pion masses, which became degenerate above the transition region, and in the 
behavior of the quark condensate which decreased to zero. 

The effective quark mass first showed a modest decrease with temperature, but above the transition region it 
started to increase and this trend was quite strong for the ZPT case. Since the condensate was already nearly zero, 
this growth was induced entirely by the meson field fluctuations. In condensed matter physics this phenomenon is 
known as pseudogap formation. At high temperatures the meson masses increased much more rapidly than the quark 
mass, so that they effectively decoupled from the system and we had a nearly ideal gas of quarks. In the ZPT case 
the quark effective mass was comparable to the temperature, however in the NOZPT case the mass was much smaller 
so that the thermodynamics closely resembled a massless quark gas. On the other hand, at low temperatures quarks 
were heavy in comparison to the temperature while pions were relatively light so that we had an ideal gas of pions 
with mass close to the physical mass. Certainly, a qualitatively similar behavior is expected on the basis of QCD. 
The transition between these two regimes occurs at a temperature which is surprisingly close to that found in lattice 
QCD simulations [T^. In addition we considered simplified approaches in which the quark mass was treated as a 
number rather than a function of the meson fields; these were found to be quite successful in the high-temperature, 
chiral-restored regime, but less accurate at low temperatures. 

We also briefly discussed the dependence of the chiral transition on the vacuum pion mass. As the chiral limit was 
approached by reducing the pion mass to a fraction < 0.75 of its true vacuum value the crossover transition became a 
first order phase transition. The critical temperature was found to be quite insensitive to the value of the pion mass, 
as in lattice analyses [T9l |. 

In the future it would be interesting to perform calculations for non-zero quark chemical potential. This would give 
the possibility of exploring the phase diagram of the model in the T — plane and comparing it with the predictions 
of other QCD-motivated models. It would also be interesting to optimize the model in order to achieve better 
correspondence with QCD. Some obvious omissions from the present approach include vector mesons, strange quarks, 
the gluon condensate and gluon-like excitations 113 ■ The latter would provide the correct degrees of freedom at 
high temperature. Also in the future, more realistic patterns of symmetry breaking, including for instance the IJ(X)a 
anomaly, should be considered. 
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APPENDIX A: CONNECTION TO PREVIOUS WORK 



The strategy followed in Refs. |21l l22l | was to assume for the purposes of counting that 

(A^) = K^> = i(A2 + 7r^)^i(^2), (Al) 

although in the numerical evaluation of the final expressions the actual values of the thermal averages were used. 
Equation (|Aip is strictly true only in the high temperature limit, where the sigma and pion masses become degenerate, 
however at low temperatures one can also show that the correct expressions are obtained. Now, writing cr^ + tt^ = 
+ 2vA, the troublesome cross term, 2vA, was expanded out. Using (|Aip the expressions needed could be cast 
in terms of a function f{v'^ -I- ■0^). 

Equation IjAip allows the general expression (|30|l to be written as a four-dimensional integral for the thermal average 
of the function in question: 

ifi-' + ^')) = J {-^) + ■ (A2) 

Since only the magnitude £^ = £o + + ^2 + ^3 occurs in the integrand, the angular integration may be carried out 
giving 

00 

{f{v'+^')) = 7^1 dU'e^v {-?pk] fiv'+e') ■ (A3) 



(■02)2 J ^ V (-0' 



Equation ljA3p then allows an easy evaluation of the principal expressions used in Refs. pTl . Defining a new 
integration variable according to = 1 -|- i'^/v'^, then setting = 2v'^/{ip'^), and integrating by parts twice one 
obtains 

1^ (^^^^) ) = 1^ (^) + 1 + (1 - ^')e''E,{z') , (A4) 

as given in 0,|23|, with Vvac denoting the vacuum value of v. The exponential integral is defined [sOl by Ei{y) — 
dtt~^e~'^* . A similar procedure yields 

iv r 



^^^^/ " ^ ~ i^')^^' erfc(z)J , (A5) 

as given in |22| . with the complementary error function defined p?3 | by erfc(z) = 1 — 2ti~^/'^ e~* dt. 

APPENDIX B: AN IDENTITY FOR THE CASE (A^) = (tt^) 

Here we show that the quark contribution to Eq. (|52(l vanishes in the case (A^) = (tt^). First we prove by induction 
that in this case 

where n and m are integers. It is simple to verify this equation in the cases m = and m = 1. After integrating by 
parts the following relations are obtained 

(A2"+2t/,2j) = (A2) [(2n + l)(A2"t/'2j) _^2j(A2"+V2j-2^j ^ 

(A2«7r202.;) = (^2^ [(A2"02,) ^ 2j {A'" TT^ ^p'^ -^)] . (B2) 
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Since (A^) = (7r|) these relations show that if Eq. IjBlj) holds for m — j — 1 then it holds for m = j and, since it holds 
for m — and 1, it therefore holds in general. The relation can also be proved by using combinatorial arguments. 
Now consider a function f{m^), where — g'^[{v + A)^ + tt^]. Expanding 2g^vA in a Taylor series 



where we have used the fact that odd powers of A give a vanishing contribution, and f^"^^ denotes the n}^ derivative 
of /. Utihzing Eq. 

(A/(m )) = 2g v2_^ ( —^^ tt^ {2g vA) 

= 2gM^V^'Hm^)) , (B4) 
again using the fact that odd powers of A give zero contribution. Thus in the case (A^) — (nf) we have the relation 



where the second equality is obtained using the relation iPHji . Equation ljB5|l was used in the text. 
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